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ABSTRACT 

We report on the analysis of a large sample of 744 type 1 Active Galactic 
Nuclei, including quasars and Seyfert 1 galaxies across the redshift range from 
< 2; < 5 and spanning nearly 6 orders of magnitude in continuum luminosity. 
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We discuss correlations of continuum and emission line properties in the rest- 
frame ultraviolet and optical spectral ranges. The well established Baldwin Effect 
is detected for almost all emission lines from O VlA1034 to [Olll] A5007. Their 
equivalent widths are significantly anti-correlated with the continuum strength, 
while they are nearly independent of redshift. This is the well known Baldwin 
Effect. Its slope P, measured as log Wx oc /3 log ALa(1450A), shows a tendency 
to become steeper towards higher luminosity. The slope of the Baldwin Effect 
also increases with the ionization energy needed to create the individual lines. 
In contrast to this general trend, the NvA1240 equivalent width is nearly inde- 
pendent of continuum luminosity and remains nearly constant. The overall line 
behaviors are consistent with softer UV continuum shapes and perhaps increasing 
gas metallicity in more luminous Active Galactic Nuclei. 

Subject headings: quasars: emission lines, galaxies: active 

1. Introduction 

Broad emission lines (BELs) are a defining property of quasar spectra. Nearly 25 years 
ago Baldwin (1977) discovered an anti-correlation between the equivalent width in C IVA1549, 
Wx{C iv), and the continuum luminosity, L;^(1450 A), measured at 1450 A in the quasar rest- 
frame. This result was based on a sample of 20 quasars spanning ~ 2 orders of magnitude 
in continuum luminosity in the redshift range 1.24<2;<3.53. It has become known as 
the "Baldwin Effect" (hereafter BEff). The BEff was subsequently confirmed in ClvA1549, 
as well as other BELs such as Lya and OviA1034 (Veron-Cetty, Veron, & Tarenghi 1983; 
Baldwin, Wampler, & Gaskell 1989; Kinney, Rivolo, & Koratkar 1990; Osmer, Porter, & 
Green 1994; Zheng, Kriss, & Davidsen 1995; Green, Forster, & Kuraszkiewicz 2001; see 
Osmer & Shields 1999 for a recent review). 

A major impetus for studying the BEff was that it might be useful for calibrating 
Active Galactic Nuclei (AGNs) luminosities, e.g. based on Wx{Cw). The AGNs could 
then be used as cosmological standard candles. But in the following years studies of bigger 
quasar samples revealed a large scatter in the anti-correlation of the continuum luminosity vs. 
emission line strength (Baldwin, Wampler, & Gaskell 1989; Zamorani ct al. 1992). The BEff 
is nonetheless important as a diagnostic of AGN structure and, perhaps, metal abundances 
(Korista, Baldwin, & Ferland 1998). The relation between the continuum luminosity and 
the relative emission line strengths and ratios can be used to study the evolution and physics 
of the quasar phenomenon (Baldwin 1999). In particular, this correlation can be used to 
test model predictions for the dependence of the shape of the continuum spectral energy 
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distribution as a function of luminosity (Binette et al. 1989; Netzer, Laor, & Gondhalekar 
1992; Zheng & Malkan 1993; Wandel 1999a,b). 

The most fundamental problem, however, is that the physical cause of the BEff remains 
unknown. It was suggested by Mushotzky & Ferland (1984) that the observed relation can 
be explained by an anti-correlation of the ionization parameter, U, and the continuum lumi- 
nosity, with U = Q{H)/47rr'^cnH, where Q(H) is the number of hydrogen ionizing photons 
emitted per second by the central continuum source, r is the distance between the contin- 
uum source and the emission line region, and uh is the hydrogen density in the line-emitting 
cloud. Assuming an additional relation of decreasing covering factor with increasing con- 
tinuum strength, the observed BEff could be well described for CivA1549. However, this 
model does not naturally explain the BEff in the full range of measured lines. It predicts, 
for example, a lack of a BEff for Lyct even though it is clearly detected (e.g., Kinney, Rivolo, 
& Koratkar 1990; Osmer, Porter, & Green 1994; Laor et al. 1995; Green 1996). 

Another important clue to the physical cause of the BEff is that the strength of the re- 
lationship (i.e., the slope of the Wx - Lc anti-correlation) seems to depend on the ionization 
energy of the emission-line species (Zheng et al. 1995; Espey & Andreadis 1999). In particu- 
lar, the equivalent widths of high-ionization lines like O VlA1034 decrease more dramatically 
with Lc than lines with moderate ionization energies like ClvA1549 or low ionization ener- 
gies like Mg IIA2798 and Balmer emission lines. Our results in the present paper confirm this 
claim, and improve upon the overall empirical characterization of the BEff correlations. 

The BEff might result from a more fundamental correlation between the continuum 

luminosity, Lc, and the shape of the ionizing (EUV-soft X-ray) continuum (Binette ct 
al. 1989; Zheng & Malkan 1993; Zheng et al. 1997; Korista ct al. 1998). Netzer (1985, 1987) 
and Netzer et al. (1992) suggested accretion disc models to explain the observed continuum 
and emission line correlations. Recently, Wandel (1999 a,b) added to this the growth in black 
hole mass by accretion in the accretion disc model. His analysis predicts that the continuum 
luminosity increases towards higher black hole mass and that the shape of the ionizing 
continuum becomes softer. Hence, it is suggested that the BEff is driven by a softening of 
the ionizing continuum towards higher luminosities. This model is attractive because the 
UV- X-ray spectral softening has been well documented by observations (Tananbaum ct 
al. 1986; Wilkes et al. 1994; Green et al. 1995), and because the spectral softening provides a 
natural explanation for steeper BEff slopes in higher ionization lines. 

An additional possibility is that the BEff is driven, at least in part, by a trend for 
higher metallicities in more luminous AGNs. Korista et al. (1998) presented theoretical 
results showing the dependence of the emission line equivalent widths on both the continuum 
shape and the metallicity of the gas. Their proposal is based on evidence for higher metal 
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abundances in more luminous quasars (Hamann & Ferland 1993, 1999; Dietrich et al. 1999; 
Dietrich & Wilhelm-Erkens 2000), and on the suggestion that more luminous quasars reside 
in more massive host galaxies, which will naturally be more metal rich (Cen & Ostriker 1999; 
Pettini 1999; Kauffmann & Haehnelt 2000; Granato et al. 2001; see Hamann & Ferland 1999 
for a review). 

Our combination of ground-based and satellite data also provides the opportunity to 
study evolutionary aspects of the line strengths in detail — e.g. by measuring the same 
rest-frame UV lines across a wide range of redshifts. One serious complication affecting 
previous work is that luminosity and rcdshift arc often correlated in quasars samples, because 
the quasars are selected from magnitude-limited surveys. Most studies have supported the 
original claim by Baldwin (1977), that scales inversely with luminosity and that there is 
no significant trend with redshift z (e.g., Kinney, Rivolo, & Koratkar 1990; Osmer et al. 1994; 
Francis & Koratkar 1995; Wilkes et al. 1999). However, some recent work based on the Large 
Bright Quasar Survey (LBQS) claims that the relationship to redshift is even stronger than 
with luminosity (Forster et al. 2001; Green, Forster, & Kuraszkiewicz 2001). 

The present paper is the first in a series in which we examine the emission-line properties 
in a large sample of 744 type 1 AGNs. The sample includes Seyfert 1 galaxies and both 

radio-loud and radio-quiet quasars spanning an unprecedented wide range in both rcdshift 
(0^-2^5) and intrinsic luminosity (~ 6 orders of magnitude). The database and data 
processing will be described in detail in a future paper. A major concern regarding BEff 
studies is that selection effects in AGN samples, e.g. pertaining to emission-line strengths, 
might bias some measurements of the BEff. Our sample has the advantage of being drawn 
from a variety of surveys, with objects selected based on radio properties, grism spectroscopy, 
or broad-band color criteria. Hence, no specific selection criteria were applied for the data. 
Another key advantage is that our sample includes a wide range of luminosities at various 
redshifts. We can therefore address the separate redshift and luminosity dependences in the 
emission line data. 

There are several interesting problems we wish to address. The most basic issue is to 
quantify the nature of empirical correlations with and/or z. In an upcoming study we 
will further examine correlations with other AGN properties such as radio-loudncss (Baldwin 
1977; Sargent, Steidel, & Boksenberg 1989; Steidel & Sargent 1991; Francis & Koratkar 
1995). We will also report on trends with or z among various emission-line metallicity 
indicators. Here we present an initial analysis on the nature of the BEff. Our approach is to 
construct composite spectra for specific Lf. and z intervals, thus providing high signal-to-noise 
spectra and allowing us to study trends in both weak and strong emission lines. 

The AGN sample is described briefly in §2 and the calculation and analysis of the 
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composite spectra in §3 and §4. The main results appear in several figures in §5. We 
compare our results with prior studies and discuss them in the context of suggested models 
in §6. Throughout this paper we use the cosmological parameters Ho = 65 kms~^Mpc^^, 
Qm = 0.3, and Q/^ = (Carroll, Press, & Turner 1992). Introducing Qx = 0.7 (Netterficld 
ct al. 2002) instead of Qx = 0.0, the luminosities at the highest rcdshifts would be ~ 10% 
smaller, while at low redshifts they would reach a maximum of ~ 20 % larger at redshift 
2; ~ 1. 



2. The Quasar Sample 

We have compiled a large sample of rcst-framc visible and ultraviolet spectra for 826 type 
1 AGNs. A majority of the spectra were obtained by several groups for different studies over 
the last 20 years using ground-based instruments as well as International Ultraviolet Explorer 
(lUE) and Huhhle Space Telescope (HST) (Bahcall et al. 1993; Baldwin et al. 1989; Chaffee et 
al. 1991; Corbin & Boroson 1996; Kinney et al. 1991; Lanzetta, Turnshek, & Sandoval 1993; 
Laor et al. 1995; Sargent et al. 1988, 1989; Schneider, Schmidt, & Gunn 1991a, b; Steidel 
1990; Steidel &: Sargent 1991; Steidel unpubhshed; Storrie-Lombardi et al. 1996; Weymann 
et al. 1991, 1998; Wills et al. 1995; Zheng et al. 1997). In total, the sample contains 351 quasar 
spectra measured with HST. Furthermore, we observed a large number of the quasars in our 
sample at redshifts -2^3 (Constantin et al. 2002; Dietrich et al. 1999; Dietrich & Wilhelm- 
Erkens 2000; Dietrich et al. 2002a, b). For this study we exclude Broad- Absorption Line 
quasars (BALQSOs), although there are indications that their emission line properties do 
not differ from non-BAL quasars (Weymann et al. 1991). We also excluded several quasar 
spectra with extremely poor signal-to-noise ratios. The sample we investigate consists of 744 
type 1 AGNs. 

Because this AGN sample was compiled from many independent projects, it encompasses 
a wide range of selection criteria. Most of the quasars especially at high redshifts were 
discovered by color selection techniques (e.g., Storrie-Lombardi et al. 1996) or found by 
objective prism surveys (e.g., Chaffee et al. 1991; Sargent et al. 1988, 1989; Schneider at 
aL 1991a, b). 

We used the radio flux densities given in Veron-Cetty & Veron (2001) to determine 
the radio loudness for the quasars. Following Kellermann et al. (1989) we calculated the 

parameter = \og{Fji/ Fb) with Fr as the radio flux density at 5 GHz and Fb as the B- 
band fiux, each in the rest frame. Our classifications of radio-loud quasars arc consistent with 
classifications available in the literature (Wills et al. 1995; Bischof & Becker 1997; Wilkes 
et al. 1999; Stern et al. 2000). In total our sample contains 265 radio-loud quasars and 479 
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radio-quiet AGNs. 

Kinney et al. (1990) and Pogge & Peterson (1992) pointed out that a significant fraction 
of the scatter in the BEff in Seyfert 1 galaxies is caused by variabihty. They found that 
observations obtained at multiple epochs for individual sources display a BEff-like correlation 
(the intrinsic BEff), but with a significantly steeper slope than the BEff observed for an 
ensemble of AGNs. The physical mechanisms which drive each effect might be different, 
i.e., the ensemble BEff might be due primarily to variations in black hole mass while the 
intrinsic BEff might be caused by changes of the accretion rate M (Green 1999). For 8 of the 
36 low luminosity AGNs (log ALa(1450A) < 44, for \L\ in crgs^-*^), the spectra arc means 
from intense monitoring campaigns (Crenshaw et al. 1996; Collier et al. 1998,2001; Dietrich 
et al. 1993; Goad et al. 1999; Kaspi et al. 1996; Korista et al. 1995; Peterson et al. 1994; 
Reichert et al. 1994; Santos-Lleo et al. 2001; Stirpe et al. 1994; Wanders et al. 1997). At the 
high luminosity end of the luminosity distribution, it has been shown that quasars vary on 
longer timescales and with smaller amplitudes (e.g., Kaspi et al. 2000). For these reasons, 
the scatter caused by intrinsic variability is not expected to strongly influence our analysis 
of the present large quasar sample, particularly when comparisons of composite spectra are 
employed. 

We transformed each quasar spectrum to the rest-frame using a redshift measured from 
ClvA1549. To determine the redshift we fit a Gaussian profile to the upper part of the 
ClvA1549 emission line having > 50 % of the peak intensity. We corrected each AGN 
spectrum for Galactic extinction (Seaton 1979), based on the neutral hydrogen column den- 
sity Nh expressed in units of 10^° cm^^, assuming Eb-v = Nh/GO (Dickey & Lockman 
1990). 

The redshift distribution of the quasar sample is presented in Figure 1 as a function 
of continuum luminosity, ALa(1450A). The quasars cover a redshift range of < 2; < 5 
and nearly 6 orders of magnitude in luminosity. For most of the redshift range a reasonable 
coverage in luminosity is achieved, spanning at least 3 orders of magnitude for ^ > 1. The 
radio- loud and radio-quiet quasars are displayed with different symbols in Figure 1. Their 
distributions across luminosity and redshift are very similar within out sample, although the 
fraction of radio-loud quasars increases for redshifts z <2. 

3. Composite Spectra 

A common approach to study the BEff is to measure the emission line properties of 
individual quasars. However, the results are often limited especially for the weaker emission 
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lines by the low signal-to-noise ratio of individual spectra. There can also be significant 
problems introduced by including upper limits on weak lines in individual spectra. To 
avoid those problems, we computed composite spectra in different intervals of redshift and 
luminosity. The analysis of these spectra has the advantage of more clearly representing 
the average properties of the sample used to calculate each composite. Figure 1 indicates 
the ranges in luminosity and redshift used to compute the composites. Except for the low 
and high luminosity ends of our sample, each composite spectrum is based on more than 20 
individual quasar spectra. The scatter introduced by different individual spectra is taken into 
account by calculating a root-mean-square (standard deviation) spectrum for each composite, 
and using that to estimate the uncertainties in the measured equivalent widths (§4). 

Before combining the individual quasar spectra, we normalized the flux densities to 
unity in a 20 A wide continuum range centered at A = 1450 A by multiplying each spectrum 
with a suitable factor. We used these normalized quasar spectra to calculate the mean and 
rms spectrum for each composite, and the number of spectra contributing in each wavelength 
bin. 

To avoid the influence of strong narrow absorption lines, we developed a search routine 
to detect strong narrow absorption features and exclude them from the calculation of the 
composite spectrum. To search for strong absorption lines we applied the following method 
(a more detailed description will be given in a subsequent paper on the properties of the entire 
quasar sample, Dietrich ct al. 2002c). First, a preliminary mean spectrum is calculated for the 
sample of spectra under study. We then divided each individual spectrum by this preliminary 
mean to derive a ratio spectrum. Each ratio spectrum is smoothed with a running boxcar 
function, which provides in addition to the mean also an rms value for each wavelength 
box. The resulting rms spectrum indicates the scatter in flux for each wavelength element 
in the immediate spectral region. A narrow absorption line will cause a sharp increase in 
the rms spectrum. The contaminated spectral regions of individual spectra, identified by 
these rms spikes, are excluded from the calculation of the final composite spectrum, and 
no interpolation has been applied. Multiple tests with difi^erent AGNs and different sample 
sizes show that this technique sucessfully removes strong narrow absorption lines from our 
composite spectra (see below) without affecting the emission lines. Significant effects of 
narrow absorption lines remain only in the Lya forest (A < 1216A) for composites that 
include quasars with redshift z > 1. 

The resulting composite quasar spectrum for the entire sample is shown in Figure 2. 
A single power law fit to the continuum at A > 1200 A is also displayed; it has a slope 
oi a — —0.43 {F^ oc v""). This spectral index is consistent with the slope reported for 
other composite spectra (Brotherton et al. 2001; VandenBerk et al. 2001; Telfer et al. 2002). 
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The bottom panel of Figure 2 shows the number of quasars contributing to each individual 
wavelength element. For every wavelength 800 < A < 5600 A this composite quasar spectrum 
is based on at least ~ 100 quasars. 

4. Emission Line Measurements 

To measure the emission line fluxes we analyzed each composite spectrum with spectral 
fitting software provided by MIDAS^. We employed a powerlaw continuum, an ultraviolet and 
optical Fe-emission template, a Balmer continuum emission template, and multiple Gaussian 
profiles for the broad emission lines in mult i- component fits using minimization. For the 
ultraviolet wavelength range we used the Fe emission template which Vestergaard & Wilkes 
(2001) carefully derived from HST spectra of IZwl (kindly provided by M. Vestergaard). 
The primary aim was to measure the broad and prominent Fe ll-emission near MgllA2798. 
However, we also used the UV template to correct for Fell and Felll emission blended 
with Clll]A1909 and other lines across the wavelength range 1400 < A < 1900 A. For the 
optical wavelength region around H/3 and H7, we used the empirical Fe ll emission template 
obtained from IZwl (kindly provided by T. Boroson; Boroson & Green 1992). We adjusted 
the line widths in these Fe emission templates to match the measured widths of G lvA1549, 
MgllA2798, or Rp. 

We calculated several bound- free Balmer continuum emission spectra for Tg = 15000 K, 
Ue = 10® — 10^° cm~^, and 0.1 < r,^ < 2 following the formalism suggested by Grandi (1982). 
These Balmer continuum emission templates were supplemented for A > 3646A with high 
order Balmer emission lines with 10 < n < 50, i.e. H?9 and higher (Storey & Hummer 1995). 
The individual high-order Balmer lines were represented by Gaussian profiles and coadded. 

The multi-component fit of a Balmer continuum, two Fe-emission templates, and a power 
law continuum yields independent measurements of the ultraviolet and optical Fe-emission 
strengths. Furthermore, the spectral index, a, of the power law continuum does not include 
contributions from Balmer continuum or from broad Fe emission. The uncertainties associ- 
ated with the fit of Balmer continuum emission and Fe emission templates were estimated 
from the distribution for each fit (cf., Dietrich et al. 2002a). 

The emission line fluxes were measured after subtracting the power law continuum, 
Balmer continuum, and Fe emission spectra. The line fluxes were determined by multi- 
component Gaussian fits. In general, we fit the emission lines with a "broad" and a "narrow" 
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Gaussian component. For C lvA1549, however, it was necessary to introduce in addition a 

very broad component (which turned out to have typically FWHM ~ 10^ kms~^). To deter- 
mine the fluxes of Lyo; and N vA1240, which are severely blended, we used the components 
obtained from the ClvA1549 profile as templates (preserving their widths and redshifts). 
Please see Figure 2 in Dietrich & Hamann (2002d) for examples of fitted profiles. Note that 
most of the Lya forest absorption on the blue side of the Lja profile is corrected by the 
methode we applied to calculated the composite spectra. The O VI A 1034 emission line is 
subject to severe corruption by Lya forest absorption. However, in composite spectra it 
can be assumed that the same fraction of the O VlA1034 emission line flux and the adjacent 
continuum flux is absorbed by the Lya forest. Hence, the equivalent width is not expected 
to be significantly altered by Lya forest absorption. The OviA1034 emission line is also 
contaminated by LjP emission in the blue wing, whose relative strength shows a wide range, 
Ly/5/0 VIA1034 = 0.30 ±0.21 (Laor et al. 1995). Our fits to the O viA1034 emission line pro- 
files were achieved with a broad and narrow Gaussian profile, and no additional component 
due to Ljf3 emission was necessary. 

The C lll]A1909 emission line complex was reconstructed with four gaussian profiles, cor- 
responding to Al IIIA1857, Si lll] A1892, and broad and narrow Gaussian profiles for C ill] A1909. 
The correction for Fe emission allows a better measurement of MgllA2798, HellA1640, 
Nlll]A1750, and the Cm] emission line complex (AllllA1857, Silll]A1892, Clll]A1909). The 
measurement of the optical Fell emission with an appropriate scaled template made it 
possible to measure the HellA4686 emission line which is otherwise severely blended with 
Fe 11(37,38) multiplets. 

The cumulative errors in the emission line flux measurements and Wx were calculated 
as the square root of the sum of i) the uncertainties in the fltting and ii) the rms scatter 
that is introduced by both the noise and distribution of the individual spectra contributing 
to that particular composite spectrum, squared. 

5. Results and Analysis 

5.1. Trends with Redshift 

To search for a dependence of emission line properties, particularly Wx, on cosmic time 
(redshift), we selected a narrow luminosity range. The location of this luminosity range 
was also selected to achieve a comparable number of individual quasars contributing to each 
redshift bin in the z - ALa(1450A) plane. To cover a redshift range as large as possible for 
an almost constant luminosity, we selected a luminosity range of 46.16 < log ALa(1450A) < 
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47.16 (i.e., 43.0 < log La(1450A) < 44.0). This range contains 323 quasars with z > 0.5 and 
an average luminosity of log ALa(1450A) = 46.67 ± 0.12. Wc computed composite spectra 
for redshift bins with Az = 0.5 centered on redshift z = 0.75, 1.25, 3.75, 4.25, and4. 75, 
respectively (Fig. I). Figures 3a and 3b show the composite spectra for each redshift bin. 
These composite spectra are normalized by the power law continuum which we derived from 
the multi-component fits in §4. 

As a first impression the relative line strengths in the normalized composite quasar 
spectra stay nearly constant with redshift z. However, some emission lines show a marginal 
trend for increasing equivalent width at the highest redshifts. To quantify this result, we 
measured the equivalent width of most emission lines in the ultraviolet wavelength range 
(AA1000-3000A) and of the broad ultraviolet Fe ll-emission feature (AA2200-3090 A). The 
results for Wx as a function of redshift are displayed in Figures 4a and 4b. In general, the 
equivalent widths of the lines remain constant within ~ 30%. The stronger emission lines, 
however, like Lya, C lvA1549, and C ill] A1909, show marginal trends for a relation of Wx and 
redshift. Wx stays nearly constant within less than ~ 20 % for redshifts less than z ~ 2 — 3, 
but then for z > 3 there are marginal indications for a slight increase of Wx at constant 
luminosity. Close inspection of Figures 4a and 4b indicates that in particular the equivalent 
widths measured for the composites with z > 4 are larger compared to lower redshifts. Most 
of the quasar spectra in our sample with z > A were obtained by Constantin et al. (2002). 
They found evidence that in spite of strong contamination by Lya-forest absorption, Lya is 
stronger for quasars with z > A than in their \ow-z counterparts. 



5.2. Trends with Luminosity 

Next we divided the quasar sample into eleven subsets with respect to the luminosity 
L;^(1450A). Except for the lowest luminosity bin, each composite spectrum covers a range 
of A log ALa = 0.5 in luminosity starting at log ALa(1450A) = 43.16 ergs~^ (Fig. 1). The 
two lowest luminosity composite spectra consist of less than 10 objects. 

Figures 5a and 5b show the composite spectra ordered by luminosity for the ultraviolet 
wavelength range (A < 2000A), while Fi gure 6 shows the longer wavelengths (AA2000 — 
5600A) . Each composite spectrum is normalized by the power law continuum fit we derived 
from the multi-component analysis (§3). These normalized spectra show very obviously a 
strong relation between the equivalent width, Wx, and the continuum luminosity represented 
by ALa(1450A), i.e. the Baldwin Effect, for most of the emission lines. Figures 7a and 7b 
show the measured values of Wx for each luminosity interval. The BEff can be seen easily in 
strong lines such as ClvA1549, O VIA1034, LyaA1216, Clll]A1909, and MgllA2798 (Figures 
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5a and 6). Taking advantage of the high signal-to-noise ratio of the composite spectra, we also 
detected this relation for weaker emission lines such as HeIlA1640. Furthermore, emission 
lines of the same ion but at different wavelengths (Olll]A1663, [Olll] A5007 and HellA1640, 
HeIlA4686) possess nearly identical slopes (Tab. 1). For the first time, we report on the 
BEff for Nlv]A1486, Olll]A1663, and Nlll]A1750. In spite of the large scatter, the broad Fe 
emission features in the ultraviolet show some indications of a BEff, while within the errors, 
the optical Fe ll emission lacks a BEff. In addition to the broad emission lines, we detect the 
Baldwin Effect also for the prominent forbidden emission fine [Olll] A5007 (Figure 6), which 
is regarded as a typical narrow-line region (NLR) emission line. In contrast to the strong 
BEff in most of the emission lines visible in Figures 5a, 5b, and 6, the Balmer lines H7 and H/3 
show only weak indications of a Baldwin Effect, which might be restricted to the line center 
of the line profiles. Close inspection of the NvA1240 emission line (Figures 5a and 8) also 
indicates that this line does not exhibit a BEff like other high ionization lines. Lines having 
the same equivalent widths in each individual spectrum will cancel out in difference spectra 
(provided that their profiles arc also the same), while lines whose equivalent widths are a 
function of luminosity will still be present in difference spectra. Note that N vA1240 has on 
average cancelled out in the difference spectra as shown in Figure 8. Hence, the difference 
spectra near Nv and the careful deblending of the LyQ;-NvA1240 emission fine complex 
both reveal that NvA1240 remains nearly constant in Wx over ~ 6 orders of magnitude in 
continuum luminosity. 



5.3. Trends with Ionization Energy 

Zheng, Fang, & Binette (1992), Espey et al. (1993), and Espey & Andreadis (1999) have 
suggested that there is a correlation between the ionization energy, Xion, a specific emission 
line and the strength (slope) of the BEff. In contrast to Espey & Andreadis (1999), wc use 
the ionization energy necessary to create the ion for collisionally excited lines like C lvA1549 
and MgllA2798. and the ionization energy needed to ionize the ion for recombination lines 
like Lya and Hell. 

The BEff slope of the lines under study shows a relatively close linear correlation with 
Xion (Figures 7a, b). We determined the slope f3, where logWA f3 ■ log ALa(1450A), of the 
BEff from a linear least square fit to logWA vs. log A La(1450A). The fits were calculated 
taking into account the full uncertainties in logW^- First, we calculated the slopes across 
the entire luminosity range, shown by the dotted line in Figures 7a and 7b. However, 
several of the emission fines indicate that the slope of the BEff is flatter at lower luminosities 
(log ALa(1450A) < 44). Therefore, we calculated additional linear flts for the luminosity 
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range log A La (1450 A) > 44 (dashed-dotted lines). It is interesting to note that the hnes 
showing a change in the BEff slope have high ionization energies with Xion ^ 25 — 30 eV. The 
slopes of the logWA vs. log A La(1450A) relations are listed in Table 1 and plotted against 
Xion in Figure 9 for the entire range (filled symbols) and for log A La(1450A) > 44 (open 
symbols). Notice, that the slope of the BEff becomes steeper for increasing Xion- However, 
as has already been shown by earlier studies, N vA1240 deviates significantly from this trend 
(Espey & Andreadis 1999). NvA1240 stays constant or shows a shghtly decreasing Wa 
towards higher luminosity. 

The slope /? for O VIA1034 is of particular interest for testing the hypothesis of a steeper 
BEff towards higher ionization energies since Xion{0~^^) is more than twice Xion{He^^) and 
nearly 15 times Xion{Mg^). Assuming a constant ratio of Ly/9/0 VlA1034 = 0.30 does not 
significantly change the slope of the BEff for OviA1034. However, /3(OviA1034) would 
become steeper by nearly 25% if Ly/3 follows the same BEff as Lya and Ly/3/LyQ;=0.059 ± 
0.04 (Laor et al 1995). 

To quantify the relationship of the BEff slopes P to Xiom we calculated hnear least square 
fits to the points in Figure 9 excluding the outlying data point for NvA1240. The resulting 
slope, rj, is the same with (77 = -0.00154 ± 0.00024) or without (77 = -0.00150 ± 0.00030) 
the OviA1034 line included. This result indicates that OviA1034 does not dominate the /3 
versus Xion relation. 

6. Discussion 

6.1. Comparisons with other Studies 

In general, the slopes of the Baldwin Effect derived here are consistent with prior inves- 
tigations. However, the large range in both luminosity and redshift that distinguishes our 
quasar sample reveals correlations which were missed in some other studies. In particular, 
we can more clearly separate the trends with L and z, and our use of composite spectra 
across a wide luminosity range more accurately probes the BEff slopes, e.g., for weaker lines. 
In the following we provide some more detailed comparisons with earlier work. 

The slopes we derived for LyaA1216, ClvA1549, Clll]A1909, and MgllA2798, are in 
good agreement with prior studies (Kinney et al. 1990; Osmer et al. 1994; Laor et al. 1995; 
Green 1996; Turnshek 1997; Zamorani et al. 1992; Zheng et al. 1995; Wang et al. 1998; Es- 
pey & Andreadis 1999). The steeper slopes for ClvA1549 in Osmer et al. (1994), Laor et 
al. (1995), and Green (1996), /3 = —0.23 compared to P — —0.13 here, might be caused by 
the luminosity range of their quasar samples, ALa(1450A) > 10^^ ergs~^. Using only the 
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W\(Clv) measurements for AL;^(1450A) > 10^^ ergs ^ in our sample, we calculate a slope 
of the Baldwin Effect of /3 = -0.20 ± 0.03 for CivA1549 (Tab. 1). 

A lack of the Baldwin Effect for ClvA1549 and Clll]A1909 is reported by Wilkes et 
al. (1999). They conclude that this is caused by the Narrow-Line-Seyfert 1 galaxies (NLSyls) 
in their sample which decrease the strength of the Baldwin Effect; a significant BEff is found 
for both Civ and Cm] if the NLSyls galaxies in their sample are omitted. We note that 
these NLSyls are also low luminosity AGN. As can be seen in Figures 5a and 5b, the Baldwin 
relation is significantly flatter for low luminosities than at higher luminosities. Any AGN 
sample favoring a narrow luminosity range can dilute or bias the derived BEff. Similarly, the 
steeper BEff slopes reported by Laor et al. (1995), for OiA1305, CiiA1335, SiivA1402, and 
HellA1640 can be attributed to their smaller luminosity range (AL;v(1350A) ~ 10^^ — 10^^ 
ergs"-*-). 

The slope for the Baldwin Effect for the high ionization emission line O VlA1034 {xion — 
114 eV) is also consistent with prior investigations. The reported slope /3(OviA1034) varies 
between -0.30 ±0.04 (Zheng et al. 1995; Turnshek 1997) and -0.15 ±0.08 (Laor et al. 1995), 
with most recent studies favoring a less steep anti-correlation (—0.18 ± 0.03; Espey & An- 
dreadis 1999). Wilkes et al. (1999) found only a marginal anti-correlation of 14/a(0 VlA1034) 
and continuum luminosity. However, they noted that this result is based on just ten ob- 
jects. Green (1996) and Green et al. (2001) did not detect a significant Baldwin Effect for 
OviA1034 in their analysis of the Large Bright QSO Survey (LBQS) data set, taking all 
quasar spectra into account, i.e. including upper limit measurements as well. However, 
using only those spectra with detected OviA1034 the Baldwin Effect is clearly present. A 
possible reason for the missing Baldwin Effect for OviA1034 in their study might be the 
very narrow luminosity range (46.2 < log AL;^(2500A) < 47.0) covered by LBQS quasars 
at these wavelengths. A comparison of our Figure 7a with Figure 2 in Green et al. (2001) 
shows that we measure the same range of OviA1034 equivalent widths for this luminosity, 
log Wa(OviA1034) ~ 1.1. Hence, we conclude that the non-detection of the Baldwin Effect 
for this high- ionization line in Green et al. (2001) is predominatly caused by their very small 
luminosity range. 

In contrast to Wills et al. (1999), we detect a significant Baldwin Effect for the promi- 
nent NLR emission line [Olll] A5007 (Fig. 7b). The slope f3 of the anti-correlation is nearly 
identical to the one which we derived for Olll]A1663 (Tab. 1, Fig. 9). Again, the lack of a 
BEff for [O III] A5007 in Wills et al. (1999) might be caused by the covered luminosity range 
of less than two orders of magnitude. 

In agreement with prior investigations, we found no significant Baldwin Effect for 
NVA1240 (Osmer et al. 1994; Laor et al. 1995; Turnshek 1997; Espey & Andreadis 1999). 
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Instead of a decreasing equivalent width for increasing continuum strength, 1^a(Nv1240) 
stays nearly constant, even though it has a high value of Xion — 78 eV. 

6.2. Comparison with Model Predictions 

Although much observational and theoretical effort has been spent to decipher the pro- 
cess which causes the Baldwin Effect, the mechanism is still unclear. Several models have 
been suggested which emphasize different aspects, i.e., contributions from optically thin 
clouds and luminosity-dependent ionization parameter, covering factor, ionizing continuum 
shape, and chemical composition of the gas. 

6.2.1. Luminosity Dependent Ionization and Covering Factor 

The dependence of the strength of the Baldwin Effect on the ionization energy indicates 
that a decreasing covering factor towards higher luminosities is not sufficient to explain 
the observed BEff for each emission line. To explain the observed BEff for ClvA1549, 
Mushotzky & Ferland (1984) presented model calculations suggesting that the BEff is caused 
by a decrease of the ionization parameter, U, as well as the covering factor towards higher 
luminosities. Although this model predicts a weaker Baldwin Effect for lower luminosity 
AGN, it also predicts that Lya will lack a BEff or show only a weak BEff. A luminosity 
dependence of the covering factor might be indicated independently by measurements of 
X-ray absorption (e.g., Lawrence & Elvis 1982). 

A variant of this picture was outlined by Shields, Ferland, & Peterson (1995), who sug- 
gested that the luminosity dependence could be caused by a luminosity-dependent covering 
factor of optically thin clouds which would emit preferentially high-ionization lines. The 
presence of this component in the BLR is suggested from several aspects of variability in 
Seyfert galaxies, as well as other arguments. While luminosity dependence of coverage by 
high-ionization clouds in the BLR is apparently still a viable scenario for understanding the 
BEff, this interpretation suffers from the lack of a strong physical basis for predicting this 
behavior; explanations are ad hoc or phenomenological at best. 

6.2.2. Luminosity dependent Spectral Energy Distribution 

Several models have been advanced that focused primarily on the spectral energy distri- 
bution (SED) of the ionizing continuum, and its consequences. This approach was motivated 



- 15 - 



by increasing evidence that the SED of the continuum becomes softer for more luminous AGN 
(e.g., Sargent & Malkan 1982; Binette ct al. 1989; Schulz 1992; Netzer et al. 1992; Zheng et 
al. 1992; Zheng & Malkan 1993; Green 1996, 1998; Wang et al. 1998). 

Empirical models have been suggested which combine a powerlaw continuum and a 
thermal UV bump. Detailed investigations of X-ray and ultraviolet continuum observations, 
i.e., the overall ionizing continuum shape, yield indications that the SED becomes softer 
for increasing luminosity. In particular, it has been shown that ttox, the two-point power 
index which connects the ultraviolet at 2500 A to the X-ray continuum at 2keV, increases 
significantly with luminosity (Tananbaum ct al. 1986; Wilkes ct al. 1994; Green et al. 1995). 
This softening can be understood by a shift of the thermal UV-bump towards longer wave- 
lengths in more luminous quasars, i.e., L^v increases more than for increasing luminosity 
resulting in a greater (cf., Binette et al. 1989; Zheng & Malkan 1993). 

Early models of geometrically thin and optically thick accretion discs focused on the 
consequences of the disc inchnation and the location of gas relative to the disc (Netzer 1985, 
1987). Although these models can not explain the BEff for several orders of magnitude in 
luminosity they provide an explanation for scatter in the Lc vs. Wx relation at a given 
luminosity. The spectral energy distribution (SED) of the accretion disc continuum depends 
on the luminosity which is in turn related to the mass of the black hole and the accretion 
rate. These dependencies were suggested as the primary mechanism to cause the Baldwin 
Effect. Recently. Wandcl (1999a,b) presented an evolutionary scenario which connects the 
growth of the black hole mass, the accretion rate and the continuum luminosity. As the 
mass of the black hole increases and consequently also the luminosity, the peak temperature 
of the UV bump decreases. This results in a shift of the peak of the thermal UV bump 
to longer wavelengths, i.e., the continuum becomes softer. Netzer et al. (1992) calculated 
the logl^;^(LyQ;) vs. logL;^(1216A) relation for a wide range of black hole masses, accretion 
rates, and inclination angles in geometrically thin accretion disc models. The spread of 
the calculated Wx(Lja) at a given luminosity is mainly caused by the different inclination 
angles of the accretion disc. Assuming that our composite spectra represent an average over 
a variety of disc inclinations, the BEff which we measured for Lya is consistent with their 
model predictions. 

The observed ionization dependence of the Baldwin Effect provides strong evidence 
for a luminosity dependent SED of the ionizing continuum and thus favors accretion disc 
models. However, the detailed models discussed by Wandel (1999b) show that the SEDs of 
accretion discs depend in a complex way on the mass of the black hole and the accretion 
scenario. In particular, the cutoff energy of the thermal UV-bump is expected to depend on 
the black hole mass. These models also predict that luminous quasars with massive black 
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holes should show low cutoff energies, while Seyfert 1 galaxies with intermediate massive 
black holes should display high cutoff energies. The low luminosity NLSyl galaxies with low 
mass black holes should also exhibit high cutoff energies. Different accretion scenarios would 
introduce additional scatter to this relation. 

In summary, accretion disc models provide a reasonable explanation of the Baldwin 
Effect, in particular for the ionization energy Xion dependent strength. Although these 
models indicate that the luminosity-dependent SED is a major mechanism which drives the 
BEff, the large possible parameter space allows for additional effects. Evidence for effects in 
addition to variations of the spectral energy distribution is given by the failure to explain 
the lack of a Baldwin Effect for another high ionization emission line, N vA1240. 

6.2.3. The Influence of Metallicity Variations 

Baldwin ct al. (1995) and Korista et al. (1997, 1998) introduced the "locally optimally- 
emitting cloud" (LOG) model, which assumes that the BLR cloud ensemble covers a wide 
range of internal densities and occurs over a wide range in distance from the central contin- 
uum source. Under these conditions the emitted spectrum is controlled by powerful selec- 
tion effects and the typical observed quasar spectrum can be naturally produced (Korista et 
al.l997, 1998). 

Korista et al. (1998) computed for a wide range of different continuum slopes {F^ oc 
i/",— 2 < a < —1), densities (rig = 10^ to lO^^cm"^), and metallicities (0.2 to 10 times 
solar) the resulting emission line spectra. In recent years there has been growing evidence 
that the chemical composition of quasar gas can reach several times solar metallicity with 
a trend of higher metaUicities in more luminous quasars (Hamann Sz Ferland 1993, 1999; 
Ferland et al. 1996, Dietrich ct al. 1999; Dietrich & Wilhclm-Erkens 2000). The emission 
lines studied by Korista ct al. (1998) follow the general trend for Wx diminishing with softer 
continuum shapes. Assuming the softer continua arc related to higher luminosities, this 
yields the Baldwin Effect. However, in contrast to those emission lines, NvA1240 lacks a 
Baldwin Effect even though its high ionization energy Xion = 78 eV suggests that the decrease 
of I^a(Nv) should be strong. Korista et al. (1998) introduced metallicity as an additional 
parameter to provide a possible explanation for the lack of a Baldwin Effect in NvA1240. 
They found in their calculations that the equivalent width of Lyct, ClvA1549, or O VlA1034 
show only a weak dependence on the gas metaUicity, but N vA1240 strongly depends on it. 

In particular, they assume that nitrogen scales like a secondary element, so that its abun- 
dance relative to the other metals increases linearly with the overall metallicity (e.g., N/0 
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oc 0/H, see Hamann et al. 2002). Korista et al. (1998) assume that the metaUicity increases 
with increasing AGN luminosity (Hamann & Ferland 1999). With this abundance behavior, 
Korista et al. (1998) showed that the equivalent width of N vA1240 is not expected to display 
a BEff. Therefore, the decreasing equivalent width that should occur in NvA1240 as part 
of the normal BEff is approximately compensated by the increasing relative abundance of 
nitrogen. The net result is that NvA1240 shows essentially no BEff. 

Overall, the BEff slopes that wc derived for the emission lines in our quasar sample are 
in good agreement with the predictions by Korista et al. (Table 1). The observed slopes of 
C IVA1549 and HcIlA1640 arc smaller than expected, while the slope in N IIl]A1750 is steeper. 
However, if only measurements with ALa(1450A) > lO^'* ergs^^ are taken into account, the 
observations yield a comparable steep anti-correlation as predicted by Korista et al. (1998). 
We will discuss these comparisons further in a forthcoming paper on the trends in metaUicity 
in this dataset (see also §6.4 below). 

6.3. Ionization Dependence and Curvature of the Baldwin Effect 

A common prediction of models that involve a softer continuum SED towards higher 
luminosities is a steeper slope of the Baldwin Effect for lines with increasing ionization 
energy Xion (Zheng et al. 1992; Zheng & Malkan 1993; Korista et al. 1998; Wandel 1999a). 
Figures 7a,b and 9 provide clear evidence for this prediction. High ionization lines like 
ClvA1549, HcllA1640, and HellA4686, and especially OviA1034 are particularly valuable. 
These emission lines respond to photoionizing energies of ~ 48 eV, ~ 54 eV, and ~ 114eV, 
respectively. Hence, they can be used to probe the location of the cutoff energy of the 
thermal UV-bump. The strong BEff in these high-ionization lines is expected if the ionizing 
continuum becomes softer with increasing continuum luminosity. In this case the relative 
number of ionizing photons with /i z/ > 50 eV decreases as the UV bump is moving towards 
longer wavelengths. While there is only a marginal trend for a shghtly steeper Baldwin 
Effect for Hell than for CiV, the higher ionization OviA1034 hne displays a significantly 
stronger BEff (Fig. 9). The decrease of the relative number of ionizing photons with a softer 
continuum, caused by the shift of the UV bump to longer wavelengths, results in weaker 
high-ionization lines compared to lower ionization ones. The shift of the UV bump to longer 
wavelengths will affect emission lines at higher ionization energy first. In addition, the 
relative continuum strength beneath the emission lines will be enhanced, starting at short 
ultraviolet wavelengths. Both effects together result in smaller W\ for higher luminosities. 

In addition to the general trend of steeper BEff for higher Xion-, the curvature or flattening 
of the BEff towards lower luminosities can be accommodated in the framework of some of the 
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luminosity-dependent SED models (e.g., Netzer et al. 1992; Wandel 1999a). The flattening 
at low luminosities has led to suggestions that Seyfert 1 nuclei do not participate in the 
Baldwin Eff'ect but have instead equivalent widths independent of luminosity (e.g., Wampler 
ct al. 1984). A close inspection of Figures 7a,b provides some evidence for the tendency of a 
flatter slope j3 of the BEff at low luminosities in comparison to higher luminosities, as can be 
seen for several emission lines (SilvA1402, Nlv]A1486, ClvA1549, HellA1640, Olll]A1663, 
Nlll]A1750, HellA4686). The turnover appears to occur near Alog La(1450A) ~ 44, in quite 
good agreement with predictions (Wandel 1999a). Within this context it is interesting to 
note that the steepening of the Baldwin Effect for log AL;^(1450A) ~ 44, is observed for 
emission lines with Xicm ^ 25 — 30 eV. 



6.4. Individual emission line pairs 

For several elements we have measured more than one emission line strength. The com- 
parison of the strength of the Baldwin Effect may provide information about the emission line 
region and the mechanism which causes the observed anti-correlation of W\ and continuum 
luminosity. 

For hydrogen we have measured two Balmer emission lines (H7, H/?) and Lya. The slope 
of the Baldwin Effect of H7 and H/3 is very similar and within the uncertainties consistent 
with zero, i.e. no Baldwin Effect for Balmer lines. The missing Baldwin Effect for these lines 
provide some argument against a luminosity dependent covering factor, at least for the gas 
component which is the dominant contributor to these emission lines. While the equivalent 
width W\{\il3) remains nearly constant, the flux ratio Ly«/H/3 decreases by a factor ~ 1.5 
from low to high luminosities (~ 8.5 to ~ 5.5). The different behaviour of the Balmer lines 
and Lyct might be also related to the complex physical processes of Lya and H/3 emission 
(Netzer et al. 1995). The prominent BEff of Lya and the lack of a BEff for H/3 may be caused 
in part by a relative stronger increase of the local continuum beneath Lya compared to the 
optical wavelength range of H/3. However, Lya shows a significant BEff which fits well in the 
slope 13 vs. Xion relation (Fig. 7a, 9). Lya might be blended with HellA1216 and Ov]A1218 
(Ferland et al. 1992; Shields et al. 1995). For standard BLR conditions the contribution of 
HeIlA1216 and Ov]A1218 are negligible (less than ~ 10% , priv.comm. K. Korista). But for 
higher densities and higher ionization parameter U, particularly for optically thin conditions, 
these emission lines can be quite strong compared to Lya (Shields et al. 1995). Because 
O v]A1218 has an ionization energy of Xion = 77.4 eV, this emission line might show a strong 
BEff. If optically thin gas is a significant fraction of the BLR gas, it might be possible 
that the moderate BEff of Lya might be caused at least in part by Ov]A1218. However, 



-19- 



if Ov]A1218 is a significant contamination of Lya it is expected from the results in Shields 
et al. (1995) that HellA1640 and ClvA1549 should be stronger than observed. In addition, 
variability behavior in some objects suggests that the optically thin gas may contribute most 
strongly to the line wings (e.g., Ferland et al. 1990; Peterson et al. 1993), but the part of the 
emission lines that changes the most with luminosity has a narrow FWHM (e.g., Osmer et 
al. 1994). Although we cannot exclude that O v]A1218 provides some contribution to Lya, 
it does appear that Lya itself shows a moderate Baldwin Effect in contrast to the Balmer 
emission lines. 

The two emission lines of He"*" which wc measured, HcIlA1640 and HcIlA4686, show 
very similar slopes for the Baldwin Effect. This behavior is expected for HciiA1640 and 
He IIA4686 as typical recombination lines. The similar slope of the BEff provides also some 
evidence that luminosity-dependent radiative transfer effects, coUisional effects or dust effects 
can not affect these emission lines. 

Another hne pair we measured is Olll]A1663 and [Olll]A5007. While the ionization 
energy necessary to create O'''^ is Xion — 35.1 eV, the excitation energies of these hues differ 
by a factor of ~ 3. If the ionizing continuum becomes softer for increasing luminosity, the 
temperature in the emission line region should drop and hence Olll]A1663 should be more 
affected than [Olll]A5007. Therefore, we might expect Olll]A1663 to have a steeper BEff 
than [OlIl]A5007. Observationally their BEffs are the same. However, the comparison of 
these two lines is very complicated because they differ in critical density by several orders of 
magnitudes (n^"* ~ 7 x 10^ cm"^ for [Olll]A5007 and n^"* ~ 3 x 10^° cm"^ for Olll]A1663). 
In particular, [Olll]A5007 forms in the narrow line region (NLR), which might be heated 
and ionized additionally by shocks, and is, in any case, spatially distinct from the BLR. 
The different kinematics of the O ill emitting gas is also given by the significantly different 
profiles we used to measure the emission line fiux and equivalent width. The Olll]A1663 
emission line was fit with a broad and narrow Gaussian profile (with an average FWHM 
of 5140 ± 460kms~^ and 1865 ± 290kms~^, respectively). The two Gaussian components 
we used to fit [OlIl]A5007 result in a profile different compared to OlIl]A1663, with average 
FWHMs of 1490 ± 50kms~i and 480 ± 25kms-^ The overall average FWHM of the fitted 
Olll]A1663 and [Olll]A5007 lines are 4130 ±450 kms^^ and 725 ± 100 kms-\ respectively 

We have studied three emission lines of nitrogen. NvA1240 lacks within the errors a 
BEff. In contrast, the intercombination lines Nlv]A1486 and NlIl]A1750 exhibit a Baldwin 
Effect that is consistent with the overall trend of f3 vs. Xion (Fig. 9). It might be possible that 
additional effects influence the emission properties of these nitrogen intercombination lines 
at higher luminosities. We will discuss these issues in our forthcoming paper on elemental 
abundances. 
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7. Conclusion 

We have investigated a large sample of 744 type 1 AGN covering the redshift range from 
< z < 5 and nearly 6 orders of magnitude in continuum luminosity. To enhance the signal- 
to-noise ratio, minimize the influence of peculiarities of individual quasars, and investigate 
weak as well as strong emission lines, we computed composite spectra representing narrow 
intervals in redshift and luminosity. The emission line fluxes were derived using multi- 
component Gaussian fits after removing a powerlaw continuum fit, a Balmer continuum 
emission template, and a UV and optical iron emission template from the composite spectra. 
Our main results are the following. 

• In composite spectra spanning the full redshift range at nearly constant luminosity we 
detect no strong trend in the line Wx with redshift i.e., with cosmic time. However, 
there is a marginal tendency for the highest redshift quasars {z > 4) to show slightly 
stronger emission lines than their counterparts at lower redshift. 

• In the composite spectra ranked by luminosity we find a significant Baldwin Effect 
in nearly all emission lines in the ultraviolet to optical domain. The only exceptions 
are NvA1240, B.(3, H7, and optical Fell which remain constant in Wx within the 
uncertanties. The lack of a BEff for the high-ionization feature NVA1240 suggests 
that the chemical composition of the gas is an additional parameter that can strongly 
influence the equivalent width of this and possibly other lines. 

• We detect a strong Baldwin Effect for the prominent NLR emission line [OlIl]A5007. 
The strength of the BEff for [Olll]A5007 is very similar to the BEff measured in 
Olll]A1663. 

• The slope /3 of the Baldwin Effect, where logl^A oc P-logXLx, shows a significant 
correlation with the ionization energy, Xion, needed to produce the lines. 

• The slope of the Baldwin Effect, /3, tends to be steeper at higher luminosities, A La(1450A) > 
lO^^ergs"^, compared to the lower luminosity regime. 

• The Baldwin Effect, its steepening towards higher luminosities, and the correlation of 
the slope (3 with Xion can all be well explained in the context of a luminosity dependent 
spectral energy distribution of the ionizing continuum. Assuming that the SED can 
be described as a combination of a powerlaw continuum and a thermal UV bump, the 
ionizing continuum becomes softer for increasing luminosity as the UV bump is shifted 
to longer wavelengths. This behaviour can be explained with accretion disc models as 
suggested by Netzer et al. (1992) and Wandel (1999a,b). 
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Fig.l — Redshift distribution of the AGN sample as a function of intrinsic luminosity 
at A = 1450 A. The open symbols represent radio-loud quasars and the filled symbols radio- 
quiet quasars. The dashed horizontal (vertical) lines indicate the luminosity (redshift) ranges 
which were used to calculate composite spectra. At the left side (top) of the figure is 
the number of the individual spectra contributing to each composite spectrum. For the 
luminosity range 46.16 < logALA(1450) < 47.16 composite spectra were calculated for the 
redshift intervals marked by the vertical dashed lines. 

Fig.2 — Composite spectrum of the entire AGN sample (top panel). The dashed line 
shows a power law fit to the continuum (1200A < A < 5800A) yielding a = —0.43 {F^ oc z/"). 
The bottom panel displays the number of AGN contributing to each wavelength element. 

Fig. 3 — (a) Normalized composite spectra are shown for redshift bins of Az = 0.5, 
starting at z = 0.5 (cf., Fig. 1) and nearly constant luminosity (46.16 < logALA(1450) < 
47.16). The spectra were normalized by a power law continuum fit. The horizontal dashed 
lines indicate the continuum level for the individual normalized composite spectra which 
were vertically shifted for better display. The normalized continuum strength is shown for 
the spectrum at the bottom of the figure. This same scale applies to all other spectra. 

(b) Same as Fig. 3a, but with an expanded vertical scale to display the dependence of 
the relative strength of weaker emission lines as a function of redshift. Strong emission fines 
with flat tops are truncated for easier display. 

Fig. 4 — (a) Line equivalent widths, Wx, vs. redshift for nearly constant intrinsic 
continuum luminosity ALa(1450A). 

(b) Same as Fig. 4a, but for additional emission lines. 

Fig. 5 — (a) Normafized composite spectra are shown for the luminosity bins (A log A La(1450A) = 
0.5 dex) as displayed in Fig. 1. The spectra were normalized with the corresponding power 
law continuum fit. The horizontal dashed lines indicate the continuum-level for the indi- 
vidual normalized composite spectra which were vertically shifted for better display. The 
normalized continuum strength is shown for the spectrum at the bottom of the figure and it 
applies for the other spectra too. 

(b) Same as Fig. 5a, but with an expanded vertical scale to display the dependence of 
the relative strength of weaker emission lines as a function of luminosity. Strong emission 
lines with fiat tops are truncated for easier display. 

Fig. 6 — Same as Fig. 5a, but the normalized composite spectra are shown at larger 
wavelengths. 
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Fig. 7 — (a) Line equivalent widths, Wx, as a function of increasing continuum luminos- 
ity AL;^(1450A). We calculated linear fits to Wx{L) for the entire luminosity range (dashed 
line), as well as a luminosities log A L;^(1450A) > 44 (dashed-dotted line). 

(b) Same as Fig. 7a, but for additional emission lines. 

Fig. 8 — Difference spectra with respect to the composite spectrum log A La(1450A) = 
47.32 ergs~^ are displayed to illustrate that PVa(N vA1240) is nearly constant, regardless how 
Wa(Nv) is measured, while the other emission lines display a prominant BEff. The dotted 
line indicates the location of NvA1240. 

Fig. 9 — The slope of the Baldwin Effect as displayed in Fig. 7a and 7b as a function 
of the ionization energy Xion needed to create the specific ions. Different lines of the same 
ion (HellA1640, HellA4686 and Olll]A1663, [Olll] A5007) show nearly identical slopes. The 
filled symbols represent the slopes based on the entire luminosity range, corresponding to 
the dashed lines in Figures 7a,b. The slopes of the Baldwin Effect for higher luminosities 
only (logALA(1450A) > 44) are plotted as open symbols. 
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Fig. 3. — (a) Normalized composite spectra are shown for redshift bins of Az = 0.5, starting 
at 2; = 0.5 (cf., Fig. 1) and nearly constant luminosity (46.1 < logALA(1450) < 47.1). The 
spectra were normalized with the corresponding power law continuum fit. The horizontal 
dashed lines indicate the continuum-level for the individual normalized composite spectra 
which were vertically shifted for better display. The normalized continuum strength is shown 
for the spectrum at the bottom of the figure and it apphes for the other spectra, too. (b) 
Same as (a), but with an expanded vertical scale to display the dependence of the relative 
strength of weaker emission lines as a function of redshift. Strong emission lines with flat 
tops are truncated for easier display. 
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Fig. 4. — Line equivalent widths, W\, vs. redshift for nearly constant intrinsic continuum 
luminosity ALa(1450A). 
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Fig. 5. — (a) Normalized composite spectra are shown for the luminosity bins 
(AlogAL;^(1450A) = 0.5 dex) as displayed in Fig. 1. The spectra were normalized with the 
corresponding power law continuum fit. The horizontal dashed lines indicate the continuum- 
level for the individual normalized composite spectra which were vertically shifted for better 
display. The normalized continuum strength is shown for the spectrum at the bottom of 
the figure and it applies for the other spectra, too. (b) Same as (a), but with an expanded 
vertical scale to display the dependence of the relative strength of weaker emission lines as a 
function of luminosity. Strong emission lines with flat tops are truncated for easier display. 
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Fig. 7. — Line equivalent widths, Wx, as a function of increasing continuum luminosity 
AL;^(1450A). We calculated linear fits to VFx(L) for the entire luminosity range (dashed 
line), as well as a luminosties log A La(1450A) > 44 (dashed-dotted hne). 
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Table 1: Slope /3 of the Baldwin Effect for the individual emission lines. A linear least 
square fit was calculated for the entire luminosity range and for log A La (1450 A) > 44; 
log Wx = a -\- 13 log ALa(1450A). In addition, for comparison theoretical predictions are 
given (Korista et al. 1998). 



line 


Xion 


entire ALA(1450A)-range 


logALA(1450A)>44 


Korista et al.(1998) 




[eV] 






/3 


OviA1034 


113.9 


-0.24 ±0.04 


-0.24 ± 0.04 


-0.22 


NVA1240 


77.7 


+0.01 ±0.01 


-0.00 ±0.02 


-0.05 


HellA1640 


54.4 


-0.16 ±0.03 


-0.20 ±0.04 


-0.17 


HellA4686 


54.4 


-0.15 ±0.04 


-0.20 ±0.09 




CivA1549 


47.9 


-0.14 ±0.02 


-0.20 ±0.03 


-0.20 


Niv]A1486 


47.4 


-0.11 ±0.04 


-0.14 ±0.08 


-0.01 


Oiii]A1663 


35.1 


-0.09 ±0.03 


-0.20 ±0.05 


-0.12 


[Oiii] A5007 


35.1 


-0.10 ±0.02 






SiivA1402 


33.5 


-0.09 ±0.01 


-0.10 ±0.02 


-0.13 


Niii]A1750 


29.6 


-0.07 ±0.03 


-0.08 ±0.04 


±0.09 


AliiiA1857 


28.4 


-0.11 ±0.03 


-0.11 ±0.03 


-0.04 


Ciii]A1909 


24.4 


-0.09 ±0.01 


-0.09 ±0.02 


-0.08 


Siiii]A1892 


16.3 


-0.16 ±0.03 


-0.16 ±0.12 


-0.06 


LyaA1216 


13.6 


-0.11 ±0.01 


-0.14 ±0.02 


-0.10 


H7A434O 


13.6 


-0.02 ±0.03 


-0.02 ± 0.03 




H/3A4861 


13.6 


-0.01 ±0.01 


±0.01 ± 0.03 




C11AI335 


11.3 


-0.10 ±0.04 


-0.01 ±0.09 


±0.00 


SiiiA1260 


8.2 


±0.02 ± 0.03 


-0.01 ± 0.09 




Fe llUV 


7.9 


-0.05 ±0.03 






Fe llopt 


7.9 


-0.02 ±0.03 






MgllA2798 


7.6 


-0.09 ±0.01 


-0.09 ±0.01 


-0.09 


O1AI305 


0.0 


-0.07 ±0.02 


-0.02 ± 0.04 


-0.05 



